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Abstract 

We present the applications of wavelet analysis meth- 
ods in constrained variational framework to calculation of 
dynamical aperture. We construct represention via exact 
nonlinear high-localized periodic eigenmodes expansions, 
which allows to control contribution to motion from each 
scale of underlying multiscale structure and consider qual- 
itative approach to the problem. 

1 INTRODUCTION 

The estimation of dynamic aperture of accelerator is real 
and long standing problem. From the formal point of view 
the aperture is a some border between two types of dynam- 
ics: relative regular and predictable motion along of ac- 
ceptable orbits or fluxes of orbits corresponding to KAM 
tori and stochastic motion with particle losses blown away 
by Arnold diffusion and/or chaotic motions. According to 
standard point of view this transition is being done by some 
analogues with map technique [1]. Consideration for aper- 
ture of n-pole Hamiltonians with kicks 
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is done by linearisation and discretization of canonical 
transformation and the result resembles (pure formally) 
standard mapping. This leads, by using Chirikov criterion 
of resonance overlapping, to evaluation of aperture via am- 
plitude of the following global harmonic representation: 
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The goal of this paper is is two-fold. In part 2 we con- 
sider some qualitative criterion which is based on more 
realistic understanding of difference between motion in 
KAM regions and stochastic regions: motion in KAM re- 
gions may be described only by regular functions (with- 
out rich internal structures) but motion in stochastic re- 
gions/layers may be described by functions with internal 
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self-similar structures, i.e. fractal type functions. Wavelet 
analysis approach [2], [3] provides us with more or less 
analytical description based on calculations of wavelet co- 
efficients/wavelet transform asymptotics. In part 3 we con- 
sider the same problem on a more quantitative level as con- 
strained variational problem and give explicit representa- 
tion for all dynamical variables as expansions in nonlinear 
periodic high-localized eigenmodes. 

2 QUALITATIVE ANALYSIS 

Fractal or chaotic image is a function (distribution), 
which has structure at all underlying scales. So, such ob- 
jects have additional nontrivial details on any level of res- 
olution. But such objects cannot be represented by smooth 
functions, because they resemble constants at small scales 
[2], [3]. So, we need to find self-similarity behaviour dur- 
ing movement to small scales for functions describing non- 
regular motion. So, if we look on a "fractal" function / 
(e.g. Weierstrass function) near an arbitrary point at differ- 
ent scales, we find the same function up to a scaling fac- 
tor. Consider the fluctuations of such function / near some 
point xo 

fioc{x) = f{x + x) - f(x ) (3) 



then we have 
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where a(xo) is the local scaling exponent or Holder expo- 
nent of the function / at Xq. 

According to [3] general functional spaces and scales of 
spaces can be characterized through wavelet coefficients 
or wavelet transforms. Let us consider continuous wavelet 
transform 

W g f(b,a)= { dx±g(—)f(x), 
JR" a \ a / 

b e R", a > 0, w.r.t. analyzing wavelet g, which is 
strictly admissible, i.e. 
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Wavelet transform has the following covariance property 
under action of underlying affine group: 



W g (\a,x Q + Xb) ~ X a{xo) W g {a,x + b) 
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So, if Holder exponent of (distribution) f(x) around the 
point x — xo is h(xo) € (n, n + 1), then we have the 



following behaviour of f(x) around x = xq: f(x) = cq + 
d(x - x ) + ...+ c n (x - x ) n + c\x - x \ h(xo) . Let 
analyzing wavelet have m (> n) vanishing moments, then 



W g (f)(x ,a) = Ca h ^W g (f)(x ,a) 
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and W g (f)(x ,a) ~ a h ^ when a -» 0. But if / e C°° 
at least in point x , then W g (/) (x a , a) <~ a™ 1 when a — > 0. 
This shows that localization of the wavelet coefficients at 
small scale is linked to local regularity. As a rule, the faster 
the wavelet coefficients decay, the more the analyzed func- 
tion is regular. So,transition from regular motion to chaotic 
one may be characterised as the changing of Holder ex- 
ponent of function, which describes motion. This gives 
criterion of appearance of fractal behaviour and may de- 
termine, at least in principle, dynamic aperture. 

3 CONSTRAINED PROBLEM FOR 
QUASI-PERIODIC ORBITS 

We consider extension of our approach [4]-[15] to the 
case of constrained quasi-periodic trajectories. The equa- 
tions of motion corresponding to Hamiltonian (1) may be 
formulated as a particular case of the general system of or- 
dinary differential equations dxi/ds = fi(xj,s), = 
1, 2n), where /» are not more than rational functions of 
dynamical variables Xj and have arbitrary dependence of 
time but with periodic boundary conditions. Let us con- 
sider this system as an operator equation for operator S, 
which satisfies the equation 
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which is polynomial/rational in x = (x\, . . . , x n , pi, . . . , 
p n ) and have arbitrary dependence on s and operator C(H, 
x, d/dt, d/dx, s), which is an operator describing some 
constraints as differential as integral on the set of dynam- 
ical variables. E.g., we may fix a part of non-destroying 
integrals of motion (e.g., energy) or areas in phase space 
(fluxes of orbits). So, we may consider our problem as 
constructing orbits described by Hamiltonian (1). In this 
way we may fix a given acceptable aperture or vice versa 
by feedback via parametrisation of orbits by coefficients of 
initial dynamical problem we may control different levels 
of aperture as a function of the parameters of the system (1) 
under consideration. As a result our variational problem is 
formulated for pair of operators (C, S) on extended set of 
dynamical variables which includes Lagrangian multipliers 
A. 

Then we use (weak) variation formulation 
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We start with hierarchical sequence of approximations 
spaces: 
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and the corresponding expansions: 
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As a result we have from (7) the following reduced system 
of algebraical equations (RSAE) on the set of unknown co- 
efficients a, of expansions (10): 



L(Sij,C k i,a,A) = 



(11) 



where operator L is algebraization of initial problem (7) 
and we need to find in general situation objects A. 
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We consider the procedure of their calculations in case of 
quasi/periodic boundary conditions in the bases of periodic 
wavelet functions with periods Tj on the interval [0,T] and 
the corresponding expansion (10) inside our variational ap- 
proach. Periodization procedure gives 



<Pj,k( x ) = "Vl<Pj,k(z ~ I) (I 3 ) 



tez 



So, tp,tp are periodic functions on the interval [0,T]. Be- 
cause tpj.k = <Pj,k' if k — /c'mod(2 J ), we may consider 
only < k < 2 J and as consequence our multiresolution 
has the form V 3 ; = L 2 [0, T] with Vj — spanl^^}^^, 1 
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[16]. Integration by parts and periodicity gives useful re- 
lations between objects (12) in particular quadratic case 

(d = d-i +d 2 ): 
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So, any 2-tuple can be represented by A^. Then our sec- 
ond (after (11)) additional algebraic (linear) problem is 
reduced according to [16] to the eigenvalue problem for 
{AJ?} < fe < 2 i by creating a system of 2- 7 homogeneous rela- 
tions in AJ? and inhomogeneous equations. So, if we have 
dilation equation in the form (p(x) — V%J2kez hkfi^x — 
k), then we have the following homogeneous relations 
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or in such form A\ d = 2 d X d , where X d = {A^} < fe < 2J . 
Inhomogeneous equations are: 
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where objects Mf{ 
cursive procedure 



< N — 2) can be computed by re- 
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So, this problem is the standard linear algebraical problem. 

Then, we may solve RSAE ([FT]) and determine unknown 
coefficients from formal expansion (10) and therefore to 
obtain the solution of our initial problem. It should be noted 
that if we consider only truncated expansion with N terms 
then we have from (|lTJ) the system of N x 2n algebraical 
equations and the degree of this algebraical system coin- 
cides with the degree of initial differential system. As a 
result we obtained the following explicit representation for 
periodic trajectories in the base of periodized (period Ti) 
wavelets (10): 



Xi {s) = XiQS) +^afy\{a), Xi{0) = Xi{Ti), (18) 
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Because affine group of translation and dilations is inside 
the approach, this method resembles the action of a mi- 
croscope. We have contribution to final result from each 
scale of resolution from the whole infinite scale of spaces. 
More exactly, the closed subspace Vj(j € Z) corresponds 
to level j of resolution, or to scale j. The solution has the 
following form 



x{s) = xf; ow (s) + J2 Zjfas), u 3 
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Figure 1: Periodic structure on level 6. 

which corresponds to the full multiresolution expansion 
in all time scales. Formula ( p^ ) gives us expansion into a 
slow part xff w and fast oscillating parts for arbitrary N. 
So, we may move from coarse scales of resolution to the 
finest one for obtaining more detailed information about 
our dynamical process. The first term in the RHS of equa- 
tion (19) corresponds on the global level of function space 
decomposition to resolution space and the second one to 
detail space. In this way we give contribution to our full 
solution from each scale of resolution or each time scale. 
On Fig. 1 we present (quasi) periodic regime on section 
x — p x corresponding to model (1). 
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